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Expressions are given for the maximum and minimum values of Poisson's ratio 
v for materials with cubic symmetry. Values less than —1 occur if and only if the 
maximum shear modulus is associated with the cube axis and is at least 25 times 
the value of the minimum shear modulus. Large values of \v\ occur in directions 
at which the Young's modulus is approximately equal to one half of its 111 value. 
Such directions, by their nature, are very close to 111. Application to data for 
cubic crystals indicates that certain Indium Thallium alloys simultaneously exhibit 
Poisson's ratio less than -1 and greater than +2. 
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1. Introduction 

The Poisson's ratio v is an important physical quantity in the mechanics of solids, 
second only in significance to the Young's modulus. It is strictly bounded between 
— 1 and 1/2 in isotropic solids, but no such simple bounds exist for anisotropic 
solids, even for those closest to isotropy in material symmetry: cubic materials. In 
fact, Ting & Chen (2005) demonstrated that arbitrarily large positive and negative 
values of Poisson's ratio could occur in solids with cubic material symmetry. The 
key requirement is that the Young's modulus in the 111— direction is very large 
(relative to other directions), and as a consequence the Poisson's ratio for stretch 
close to but not coincident with the 111— direction can be large, positive or negative. 
Ting & Chen's result replaces conventional wisdom (eg. Baughman et al. 1998) 
that the extreme values of v are associated with stretch along the face diagonal 
(110— direction). Boulanger & Hayes (1998) showed that arbitrarily large values 
of |^| are possible in materials of orthorhombic symmetry. Both pairs of authors 
analytically constructed sets of elastic moduli which show the unusual properties 
while still physically admissible. The dependence of the large values of Poisson's 
ratio on elastic moduli and the related scalings of strain are discussed by Ting 
(2004) for cubic and more anisotropic materials. 

To date there is no anisotropic elastic symmetry for which there are analytic 
expressions of the extreme values of Poisson's ratio for all materials in the sym- 
metry class, although bounds may be obtained for some specific pairs of directions 
for certain material symmetries. For instance, Lempriere (1968) considered Poi- 
son's ratios for stretch and transverse strain along the principal directions, and 
showed that the is bounded by the square root of the ratio of principal Young's 
moduli, |f(n, m)| < (E(n)/E(m)) 1 / 2 (in the notation defined below). Gunton and 
Saunders (1975) performed some numerical searches for the extreme values of v in 
materials of cubic symmetry. However, the larger question of what limits on v exist 
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for all possible pairs of directions remains open, in general. This paper provides 
an answer for materials of cubic symmetry. Explicit formulas are obtained for the 
minimum and maximum values of v which allow us to examine the occurrence of 
the unusually large values of Poisson's ratio, and the conditions under which they 
appear. Conversely, we can also define the range of material parameters for which 
the extreme values are of "standard" form, i.e. associated with principal pairs of 
directions such as u(110, 110) for stretch and measurement along the two face di- 
agonals. For instance, we will see that a necessary condition that one or more of 
the extreme values of Poisson's ratio is not associated with a principal direction is 
that f (110, 110) must be less than —1/2. The general results are also illustrated by 
application to a wide variety of cubic materials, and it will be shown that values of 
v < — 1 and v > 2 are possible for certain stretch directions in existing solids. 

We begin in §2 with definitions of moduli and some preliminary results. An 
important identity is presented which enables us to obtain the extreme values of 
both the shear modulus and Poisson's ratio for a given choice of the extcnsional 
direction. Section 3 considers the central problem of obtaining extreme values of 
v for all possible pairs of orthogonal directions. The solution requires several new 
quantities, such as the values of v associated with principal direction pairs. Section 
4 describes the range of possible elastic parameters consistent with positive definite 
strain energy. The explicit formulae the global cxtrcma are presented and their 
overall properties are discussed in §5. It is shown that certain Indium Thallium 
alloys simultaneously display values of v below —1 and above +2. 

2. Definitions and preliminary results 

The fourth order tensors of compliance and stiffness for a cubic material, S and 
C = S _1 , may be written (Walpole 1984) in terms of three moduli k, and yu 2 , 

S* 1 - (3k)* 1 J + (2 Ml )^ (I - D) + (2 A * 2 ) =F1 (D - J). (2.1) 

Here I ijk i = ^(SikS^i + SuSjk) is the fourth order identity, J ijkl = \5ij5u, and 

Dijki = Su6ji§kiSii + S i2 5j2Sk2Si2 + 8i 3 Sj 3 Sk3Si 3 . (2.2) 

The isotropic tensor J and the tensors of cubic symmetry (I — D) and (D — J) arc 
positive definite (Walpole 1984), so the requirement of positive strain energy is that 
k, ii\ and fi2 are positive. These three parameters, called the "principal elasticities" 
by Kelvin (Thomson 1856), can be related to the standard Voigt stiffness notation: 
k = (en + 2ci2)/3, hi = C44 and [i2 = (en — Ci 2 )/2. Alternatively, k = (sn + 
2si2) _1 /3, fii — s 44 and fj-2 — (sn — si2) _1 /2 in terms of the compliance. 

Vectors, which are usually unit vectors, are denoted by lowercase boldface, e.g. 
n. The triad {n, m, t} represents an arbitrary orthonormal set of vectors. Directions 
are also described using crystallographic notation, e.g. n = 110 is the unit vector 
(l/\/2, — 1/V2, 0). The summation convention on repeated indices is assumed. 

(a) Engineering moduli 

The Young's modulus E(n) sometimes written E n , shear modulus G(n, m) and 
Poisson's ratio v(n, m) are (Hayes 1972) 

E(n) = l/s' u , G(n,m) = l/4, u(n, m) = -s' 12 /s' u , (2.3) 
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where s' n = s ijk ininjn k ni, s' 12 = s^kin^jinkmi and s' 44 = A ts i]k i'nim 2 n k mi. Thus, 
E(n) and v(n, m) are defined by the axial and orthogonal strains in the n and m 
directions, respectively, for an applied uniaxial stress in the n direction. E and G 
are positive while v can be of either sign or zero. A material for which v < is called 
auxetic, a term apparently introduced by K. Evans in 1991. Gunton and Saunders 
(1975) provide an earlier but informative historical perspective on Poisson's ratio. 
Love (1944) reported a Poisson's ratio of "nearly —1/7" in Pyrite, a cubic crystalline 
material. 

The tensors I and J are isotropic, and consequently the directional dependence 
of the engineering quantities is through D. Thus, 

- - — + — -(—-— )F(n), (2.4) 
E 9k 3^2 >2 Mi 

4 = - + (---)2£>(n,m), (2.5) 

- = -— + — -(—-— )-D(n,m), (2.6) 
E 9k 6[i 2 >2 IH } 2 V ' h V 1 

where 

TTi/ \ 22,22,22 rM \ 22,22,22 In n\ 

b (n) = n 1 n 2 + n 2 n 3 + n 3 n x , E>(n, m) = n 1 m 1 + n 2 m 2 + n 3 m 3 . (2.7) 
We note for future reference the relations 

D(n, m) + £>(n, t) = 1 - {n\ + n% + n\) = 2F(n) . (2.8) 



(b) General properties of E, G and related moduli 

Although interested primarily in the Poisson's ratio, we first discuss some general 
results for E, G and related quantities in cubic materials: the area modulus A, 
and the traction-associated bulk modulus K, defined below. The extreme values 
of E and G follow from the fact that < F < 1/3 and < D < 1/2 (Walpole 
1986; Hayes & Shuvalov 1998). Thus, G min ,ma X = min, max(/ii, fj, 2 ), E mi ^ maK = 
3 [(3KT 1 +G m ; n max ] , and E min ,E mlix = E 00 i,E n i for ^ > /i 2 , with the values 
reversed for < [i 2 (Hayes & Shuvalov 1998). As noted by Hayes & Shuvalov 
(1998), the difference in extreme values of E and G are related by 

3/ E m [ n 3/i? max = 1/G m i n 1/G max . (2-9) 

The extreme values also satisfy 

^/-^mir^max ^/G m m,max 1/(3k). (2.10) 

The shear modulus G achieves both minimum and maximum values if n is directed 
along face diagonals, that is, G m i n < G < G max for n = 110. 

The area modulus of elasticity A(n) for the plane orthogonal to n is the ratio 
of an equibiaxial stress to the relative area change in the plane in which the stress 
acts (Scott 2000). Thus, 1/A(n) = Sijki(Sij — niUj){8ki — n k ni) . Using the equations 
above it may be shown that, for a cubic material, 

1/A(n) - 1/E(n) = 1/(3k) . (2.11) 
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The averaged Poisson's ratio v(n) is defined as the average over m in the orthogonal 
plane, or z7(n) = [^(n, m)+z/(n, t)]/2. The following result, apparently first obtained 
by Sirotin & Shaskol'skaya (1982), follows from the relations (2.8), 

[1 - 2F(n)] /E(n) = 1/(3k) . (2.12) 

Equation (2.12) indicates that the extrema of F(n) and E(n) coincide. The traction- 
associated bulk modulus K(ii), introduced by He (2004), relates the uniaxial stress 
in the n direction to the relative change in volume in anisotropic materials. It is 
defined by 3K(n) — 1/ 'sakinkni, and for cubic materials is simply K(n) = k. It is 
interesting to note that the relations (2.10) through (2.12) have the same form as 
for isotropic materials, for which E, G, v, A and K are constants. Equations (2.4) 
to (2.6) imply other identities, e.g. that the combination 1/G + Av/E is constant. 

Further discussion of the extremal properties of G and v requires knowledge 
of how they vary with m for given n, and in particular, the extreme values as a 
function of m for arbitrary n, considered in the next subsection. Note that n = 1 1 1 
and n = 001 are the only directions for which i/(n, m) and G(n, m) are independent 
of m. It will become evident that n = 111 is a critical direction, and we therefore 
rewrite E and v in forms emphasizing this direction: 

1 - 1 '[k-mx. ^ = ^ + [i- D (n,m)]i ,2.13) 



E(n) Em L 3 y u A ' E(n) E lu L 3 v ' /J 2 

where £ m = £7(111), v\\\ = ^(111, •) and \ (Hayes & Shuvalov 1998) are 

E U1 =(— + —) , v ul = - —f- , X = • 2.14 

Both Em and v\\\ are independent of //2- The fact that F < 1/3 with equality for 
n = 111 implies that this is is the only stretch direction for which E, and hence v, 
are independent of /i 2 . Equations (2.13) indicate that E(n) and i/(n, m) depend on 
fi 2 at any point in the neighbourhood of 111, with particularly strong dependence 
if Hi is small. This singular behaviour is the reason for the extraordinary values of 
v discovered by Ting & Chen (2005) and will be discussed at further length below 
after we have determined the global extrema for v. 

(c) Extreme values of G and v for fixed n 
For a given n consider the defined vector 

m(A)E/, (^A'^'^' (2 ' 15) 

with p chosen to make m a unit vector. Requiring n • m = implies that m(A) is 
orthogonal to n if 

' !i "* 0, (2.16) 



n\ — A n\ — A n\ — A 
i.e., if A is a root of the quadratic 

A 2 - 2\{n\n\ + n\n\ + n\n\) + Zn\n\n\ = 0. (2.17) 

It is shown in Appendix A that the extreme values of D(n, m) for fixed n coincide 
with these roots, which are non-negative, and that the corresponding unit m vectors 
provide the extremal lateral directions. The basic result is described next. 
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A fundamental result: 

Let < A_ < A + < 1/2 be the roots of (2.17) and m_, m + the associated 
vectors from (2.15), i.e. 



A± ={n\n\ + -n\n\ + n\n\) ± J {n\n\ + n\n\ + nln\) 2 - Zn\n\n\, (2.18a) 



m± =p± (— , 2 r - , — r— J, 2.18b) 

r n i | "1 | n l 1-1/2 

P±=[ R^F + R^F + R^F ] • (2 ' 18c) 

The extreme values of D for given n are A± associated with the orthonormal triad 
{n, m_, m + }, i.e. 

D min (n) =D(n,m_) = A_, £> max (n) = D(n, m+) = A+. (2.19) 

The extreme values of G and ^ for fixed n follow from cqs. (2.5) and (2.6). 

The above result also implies that the extent of the variation of the shear mod- 
ulus and the Poisson's ratio for a given stretch direction n are 

1/G min (n) - 1/G max (n) = | X | 4H(n), (2.20a) 
*w(n) - f mi „(n) = |xl E(n)H(n) , (2.20b) 

where H (n) is, see figure 1, 

H(n) = [(n\n\ + n\n\ + njnl) 2 - 3n?n|n^] 1/2 . (2.21) 



3. Poisson's ratio 

We now consider the global extrema of t/(n, m) over all directions n and m. Two 
methods are used to derive the main results. The first uses general equations for a 
stationary value of v in anisotropic media to obtain a single equation which must 
be satisfied if the stationary value lies in the interior of the triangle in figure 2. 
It is shown that this condition, which is independent of material parameters, is 
not satisfied, and hence all stationary values of v in cubic materials lie on the 
edges of the triangle. This simplifies the problem considerably, and permits us to 
deduce explicit relations for the stationary values. The second method, described 
in Appendix B, confirms the first approach by a comprehensive numerical test of 
all possible material parameters. 

(a) General conditions for stationary Poisson's ratio 

General conditions can be derived which must be satisfied in order that Poisson's 
ratio is stationary in anisotropic elastic materials (Norris 2006). These are: 

s' 14 = 0, 2us' lb + s' 25 = 0, (2i/ - l)«i 6 + 4b = 0, (3.1) 

where the stretch is in the 1' direction (n) and 2' is the lateral direction (m). 
The conditions may be obtained by considering the derivative of v with respect to 
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Figure 2. The irreducible l/48th of the cube surface is defined by the isosceles triangle 
with edges 1, 2 and 3. The vertices opposite these edges correspond to, n = 111, 110 and 
001, respectively. Note that the edge 3' is equivalent to 3 (which is used in Appendix B). 



rotation of the pair (n, m) about an arbitrary axis. Setting the derivatives to zero 
yields the stationary conditions (3.1). 

The only non-zero contributions to s' 14 , s' 15 , s' 25 , s' 16 and s 26 in a material of 
cubic symmetry come from D. Thus, we may rewrite the conditions for stationary 
values of v in terms of D' u = D' 1123 etc., as 

D' u = 0, 2vD' 15 +D' 25 = 0, (2p-1)D' 16 + D' 26 =0. (3.2) 

The first is automatically satisfied by virtue of the choice of the direction m as 
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either of m± . Regardless of which is chosen, 



D' 14 = Dijkiriinjm+km- 



= P+P- 
= 0. 



+ 



(n?-A+)(n?-A_) {n\ - A+)(n| - A_) (n§ - A + )(n§ - A_) 



(3.3) 



The final identity may be derived by first splitting each term into partial fractions 
and using the following (see Appendix A) 



1. 



(3.4) 



A± — A± rig — A± 

With no loss in generality consider the specific case of m = m a , A = A a where 
a = + or a = — , and in either case, b = —a. It may be shown without much 
difficulty (Appendix B) that p± > for n in the interior of the triangle of figure 2. 
It then follows that inside the triangle, 



D' D\ 



__L 5 

Pb 



16 

Pa 



= 1, 



Pb 



A a - Ab 



^6 
Pa 



A a - A 6 



2. 



(3.5) 



These identities may be obtained using partial fraction identities similar to those 
in eqs. (3.3) and (3.4). Equations (3.2) 2 and (3.2) 3 can be rewritten 



^15 



D> 



^16 u 26 



25 

D 



16 



2v 
1 



However, using (3.5), the determinant of the matrix is 

D[ 5 D' 26 - (D[ 5 + D' 25 )D' W = -3p+p-, 



(3.6) 



(3.7) 



which is non-zero inside the triangle of figure 2. This gives us the important result: 
there are no stationary values of v inside the triangle of figure 2. Hence, the only 
possible stationary values are on the edges. 



(b) Stationary conditions on the triangle edges 

The analysis above for the three conditions (3.2) is not valid on the triangle 
edges in figure 2 because the quantities p± become zero and careful limits must be 
taken. We avoid this route by considering the conditions (3.2) afresh for n directed 
along the three edges. We find, as before, that D' 14 = on the three edges, so that 
(3.2)! always holds. Of the remaining two conditions, one is always satisfied, and 
imposing the other condition gives the answer sought. 

The direction n can be parameterized along each edge with a single variable. 
Thus, n = lpO, < p < 1, on edge 1. Similarly, edges 2 and 3 are together covered 
by n = lip, with < p < oo. In each case we also need to consider the two possible 
values of m, which we proceed to do, focusing on the conditions (3. 2)2 and (3.2)3. 
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(i) Edge 1: n = lpO, < p < 1, and m = plO or 001 

For m = plO we find that D' 15 = D' 25 = 0, D[ 6 = —D' 2(> = p—p 3 . Hence equation 
(3.2) 2 is automatically satisfied, while equation (3.2) 3 becomes 

(v-l)(p-p 3 ) = 0. (3.8) 

Conversely, for m = 001 it turns out that D' le = D' 26 = 0, and D' 15 = —D' 25 = p—p 3 . 
In this case the only non-trivial equation from equations (3.2) is the second one, 

v(p-p 3 )=0. (3.9) 

Apart from the specific cases v = or v — 1, equations (3.8) and (3.9) imply that 
stationary values of v occur only at the end points p = and p = 1. Thus, i^(001), 
v(110, 110), and ^(110,001) are potential candidates for global extrema of v. 

(ii) Edges 2 and 3: n = lip, < p < 00 and m = 110 

Proceeding as before we find that D' w = D' 26 = 0, D' 15 = \f2p(\ -p 2 )/(2+p 2 ) 2 , 
and D' 25 = p/[V2(2 + p 2 )]. Hence equation (3.2)3 is automatically satisfied, while 
equation (3.2) 2 becomes 

p[(l - Av)p 2 + 2 + Av] = . (3.10) 

The zero p — corresponds to n = 110 which was considered above. Thus, all three 
conditions (3.2) are met if p is such that 

p 2 = (^ + l/2)/(^-l/4). (3.11) 

Further progress is made using the representation of equation (2.13) combined 
with the limiting values of D which can be easily evaluated. We find 

„ (11Pll I„)|!I = mi _l ( I^ )W (3.12b) 

Substituting for p 2 from equation (3.11) into (3.12) gives two coupled equations for 
E llp and j/(llp,110): 

1 1 , x < -in , x (313) 



Eu p Em A&v 2 ' E llp Em 24v ' 
Eliminating En p yields a single equation for possible stationary values of ^(llp , 110): 

v 2 -vvm- ^-EmX = 0- (3.14) 
We will return to this after considering the other possible m vector. 
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(iii) Edges 2 and 3: m_ = pp2 

In this case D[ 5 = D' 25 = 0, D' w = y/2p(\ - p 2 )/(2 + p 2 ) 2 , and D' 26 = p(p 2 - 
4)/[\/2(2+p 2 ) 2 ]. Equation (3.2) 2 holds, while equation (3.2) 3 is zero if p = 0, which 
is disregarded, or if p is such that 

p 2 = (i/-3/2)/(i/-3/4). (3.15) 

The Young's modulus is independent of m and given by (3.12a), while v satisfies 

fit ^^h^ (1 -p 2 )(4-p 2 ) 

" (UP ' PP2) E^ p = 1/111 + 6(2 +p 2 ) 2 ElllX ■ (3 - 16) 

Using the value of p 2 from (3.15) in equations (3.12a) and (3.16) yields another pair 
of coupled equations, for E\ lp and z^(llp,001): 

i i , x » , xiy-\) (317) 



E Up E U1 48(is - l) 2 ' E Up E in 24(v - I) 2 ' 
These imply a single equation for possible stationary values of f (lip , 001): 

{y - l) 2 - {y - l)(is ln - 1) - ^niX = . (3.18) 

(c) Definition of v\ and v 2 

The analysis for the three edges gives a total of seven candidates for global ex- 
trema: i/(001), z/(110, 110), and ^(110,001) from the endpoints of edge 1, and the 
four roots of equations (3.14) and (3.18) along edges 2 and 3. The latter are very 
interesting because they are the only instances of possible extreme values associated 
with directions other than the principal directions of the cube (axes, face diago- 
nals) . Results below will show that five of the seven candidates are global extrema, 
depending on the material properties. These are ^(001), ^(110, 110), ^(110,001) 
and the following two distinct roots of equations (3.14) and (3.18), respectively, 



-i = ^m-^fu + ^m + lH^-l), (3.19a) 



V2 = ^m + l) + W(-in-l) 2 + i^m + l)(— ( 3 ' 19b ) 

The quantity £dnX nas been replaced to emphasize the dependence upon the two 
parameters v X \\ and the anisotropy ratio fJ.i/^2- The associated directions follow 
from equations (3.11) and (3.15), 

i/j = i/(ll Pl , 110), pr = ( l^yl ^j ' > (3.20a) 

V2 = V(UP2,P2P22), P2 = ( ^3/4 ) • ( 3 - 20b ) 

A complete analysis is provided in Appendix B. At this stage we note that v\ is 
identical to the minimum value of v deduced by Ting & Chen (2005), i.e. equations 
(4.13) and (4.15) of their paper, with the minus sign taken in equation (4.13). 
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4. Material properties in terms of Poisson's ratios 

Results for the global extrema are presented after we introduce several quantities. 

(a) Nondimensional parameters 

It helps to characterize the Poisson's ratio in terms of two nondimensional ma- 
terial parameters which we select as u and xo, where 

vq = -si 2 /*ii = (3/s-2/i 2 )/(6K + 2/i2) ) (4.1) 
Xo = (2 S n-2 Sl2 - S44 )/ Sn = ( ^; ( ^ ) _ 1 . (4.2) 

That is, vq is the axial Poisson's ratio z/(001, •), independent of the orthogonal 
direction, and xo = x/ s n is the nondimensional analogue of \- Thus, 

Kn ' m) = l-xoF(n) ' ^ 

a form which shows clearly that ^ is negative (positive) for all directions if v$ < 
and xo > {vq > and % < 0). These conditions for cubic materials to be 
completely auxetic (non-auxetic) were previously derived by Ting & Barnett (2005). 
The extreme values of the Poisson's ratio for a given n are 

where F is defined in (2.7) and H in (2.21). Thus, v + is the minimum (maximum) 
and V- the maximum (minimum) if xo > ( Xo < 0), respectively. 

The Poisson's ratio is a function of the direction pair (n, m) and the material 
parameter pair (^o,Xo)> i- e - v — v(n, m, vq, xo). The dependence upon v has an 
interesting property: for any orthonormal triad, 

i/(n, m, v , Xo) + v{n, t, 1 - u , xo) = 1 • (4.5) 

This follows from (4.3) and the identities (2.8). Result (4.5) will prove useful later. 
Several particular values of Poisson's ratio have been introduced: vq = ^(001, m), 
= v (111, m) associated with the two directions 001 and 111 for which v is inde- 
pendent of m. These are two vertices of the triangle in figure 2. At the third vertex 
(n = 110 along the face diagonals) we have ^(110, m) = mj^ooi + (1 — "^D^iio 
where, in the notation of (Milstein & Huang 1979), i^ooi = ^(110, 001) and = 
i/(110, 110). Three of these four values of Poisson's ratio associated with princi- 
pal directions can be global extrema, and the fourth, i/ m plays a central role in 
the definition of v\ and V2 of (3.19). We therefore consider them in terms of the 
nondimensional parameters v and xo : 



- gXo v vq - 4X0 , . 

^111 = -; r — ' ^001 = 1 — , viio = - — r — ■ ( 4 - 6 ) 

I-5X0 I-3X0 1 - 4X0 



We return to v\ and V2 later. 
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1.5r 




-2 1 1 1 1 1 1 1 ' 

-1.5 -1 -0.5 V 0.5 1 1.5 2 2.5 

001 

Figure 3. The interior of the triangle in the vooi, Vila plane represents the entirety of 
possible cubic materials with positive definite strain energy. The vertices correspond to 
k = 0, /J,i = and /12 = 0, as indicated. The edges of the triangle opposite the vertices 
are the limiting cases in which kT 1 , /i^ -1 and jx^ 1 vanish, respectively. The dashed curves 
correspond to uo — (vertical) and xo — (diagonal) and the regions a, b, c and d defined 
by these lines coincide with the four cases in equation (4.9) respectively. 



(&) Positive definiteness and Poisson's ratios 

In order to summarize the global extrema on v we first need to consider the 
range of possible material parameters. It may be shown that the requirements for 
the strain energy to be positive definite: n > 0, /i 2 > and fix > 0, can be expressed 
in terms of v§ and xo as 

-1<i*<1/2, X o< 2(l + i/o). (4.7) 

It will become evident that the global extrema for v depend most simply on the 
two values for n along a face diagonal: 1/001 & n d V\\q ■ The constraints (4.7) become 

-l<z/ 1I o<l, --(1 - I/no) < fooi < 1 - v\ia , (4.8) 

which define the interior of a triangle in the 1/001 , ^110 plane, see figure 3. This figure 
also indicates the lines Vq — and \o — (isotropy). It may be checked that the 
four quantities {i>q, vm, vqoi, ^110} are different as long xo 7^ 0> with the exception 
of vqoi and vq which are distinct if ^0X0 7^ 0. Consideration of the four possibilities 
yields the ordering 

< nooi < ^0 < ^111 < "no < 1 i0r ^0 > 0, xo < 0, (4.9a) 

-1 < i/ < ^001 < fni < "no < for i/ < 0, xo < 0, (4.9b) 

-1 < i/ ll0 < fm < vooi < < for u < 0, xo > 0, (4.9c) 

< v m < vm < u < vooi < 2 for u > 0, xo > 0. (4.9d) 
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Table 1. The global minimum of Poisson's ratio for cubic materials 





n 


m 


condition 1 


condition 2 


Fig. 4 


< vooi 


110 


001 


vooi > 


"110 > "001 


a 


-\ < "llO 


110 


110 


"llO > ~\ 


"110 < "001 


b 


-1 < "o 


001 


arbitrary 


"ooi < 


"110 > "001 


c 


— OO < V\ 


n Pl 


110 


"110 < ~\ 


"110 < "001 


d 




Table 2. 


The global 


maximum of Poisson's ratio 




^raax 


n 


m 


condition 1 


condition 2 


Fig. 5 


V 1 < -i 


llPx 


Ho 


"110 < ~\ 


"110 > "001 


a 


v < 


001 


arbitrary 


"ooi < 


"110 < "001 


b 


"llO < 1 


110 


Ho 


"110 > ~\ 


"110 > "001 


c 


VQOl < | 


110 


001 


< z/ooi < | 


"110 < "001 


d 


V2 < OO 


llp 2 


P 2 P 2 2 


"001 > | 




e 



Note that vm is never a maximum or minimum. We will see below that (4.9a) is the 
only case for which the extreme values coincide with the global extrema for v. This 
is one of the reasons the classification of the extrema for v is relatively complicated, 
requiring that we identify several distinct values. In particular, the global extrema 
depend upon more than sgn^o and sgnxo, but are best characterized by the two 
independent nondimensional parameters fooi and vn . 
We are now ready to define the global extrema. 

5. Minimum and maximum Poisson's ratio 

Tables 1 and 2 list the values of the global minimum j/ m j n and the global maximum 
Vmax, respectively, for all possible combinations of elastic parameters. For table 1, 
z/(001), z/qoi and v x \q arc defined in (4.6), and v\ and p\ are defined in (3.19a) 
and (3.20a). For table 2, v 2 and p 2 are defined in (3.19b) and (3.20b). No second 
condition is necessary to define the region for case e, which is clear from figure 5. 
The data in tables 1 and 2 are illustrated in figures 4 and 5, respectively, which 
define the global extrema for every point in the interior of the triangle defined by 
(4.8). The details of the analysis and related numerical tests leading to these results 
are presented in Appendix B. 

(a) Discussion 

Conventional wisdom prior to Ting & Chen (2005) was that the extreme values 
were characterized by the face diagonal values ^ooi and vn . But as equation (4.9) 
indicates, even these are not always extrema, since vq ~ j/(001,m) can be maxi- 
mum or minimum under appropriate circumstances (equations (4.9c) and (4.9d), 
respectively). The extreme values in equation (4.9) are all bounded by the limits of 
the triangle in figure 3. Specifically, they limit the Poisson's ratio to lie between —1 
and 2. Ting & Chen (2005) showed by explicit demonstration that this is not the 
case, and that values less than —1 and larger than 2 are feasible, and remarkably, 
no lower or upper limits exist for v. 
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-1.5 -1 -0.5 v 0.5 1 1.5 2 2.5 

001 



Figure 4. The global minimum of Poisson's ratio based on table 1. The value of v m i n 
depends upon the location of the cubic material parameters in the four distinct regions 
a, b, c and d, defined by the heavy lines inside the triangle of possible materials. The 
diagonal dashed line delineates the region in which Vmi n < —1, from eq. (5.1). 

The Ting & Chen "effect" occurs in figure 4 in the region where v m \ n = v\ and 
in figure 5 in the region ^ max = ^2- Using equation (3.19a) we can determine that 
v m\vL is strictly less than —1 if (fJ-i/fJ-2 — 1) > 24. Similarly, equation (3.19b) implies 
that v max is strictly greater than 2 if (/ii//i2 — l)(^iii + 1) > 24(2 — i^m). By 
converting these inequalities we deduce 

fmin < -1 •£> fJ.2 < 777 ^ ^001 - 13^H > 12, (5.1a) 

25 

i> max >2 <S> V2<(— + —y 1 & 13MB01 - I'lIO > 24. (5.1b) 

Ml K 

The two subregions defined by the vqoi, Viiq inequalities are depicted in figures 4 
and 5. They define neighbourhoods of the ^2=0 vertex, i.e. (I'ooij ^no) = (2. — 1), 
where the extreme values of v can achieve arbitrarily large positive and negative 
values. The condition for v mm < —1 is independent of the bulk modulus k. Thus, 
the occurrence of negative values of v less than — 1 does not necessarily imply that 
relatively large positive values (greater than 2) also occur, but the converse is true. 
This is simply a consequence of the fact that the dashed region near the tip /12 = 
in figure 5 is contained entirely within the dashed region of figure 4. 

These results indicate that the necessary and sufficient condition for the oc- 
currence of large extrema for v is that /i2 is much less than either fii or n. fi 2 is 
either the maximum or minimum of G, and it is associated with directions pairs 
along orthogonal face diagonals, ^2 = G(110, 110). Hence, the Ting & Chen effect 
requires that this shear modulus is much less than \i\ = G(001, m), and much 
less than the bulk modulus k. In the limit of very small [12 equations (3.19) give 
v\.2 ~ TvO^ni + l)/ i i/(24/i2)- Ting (2004) found that the extreme values are 
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v w ±y / 3/(165)+0(l) for small values of their parameter (5. In current notation 
this is 5 = 9/[l + % xL an d replacing Bm x the two theories are seen to agree. 

The implications of small fi2 for Young's modulus are apparent. Thus, E m - U1 / E max 
=0(^2/^1), and equation (2.13)i indicates that E(n) is small everywhere except 
near the 111 direction, at which it reaches a sharply peaked maximum. Cazzani & 
Rovati (2003) provide numerical examples illustrating the directional variation of 
E for a range of auxetic and non-auxetic cubic materials, some of which are consid- 
ered below. Their 3-dimcnsional plots of E(n) for materials with very large values 
of (ii/fjQ (see Table 3 below) look like very sharp starfish. Although the directions 
at which v\ and v^. are large in magnitude are close to the 111 direction, the value of 
E in the stationary directions can be quite different from Em. The precise values 
of the Young's modulus, Eu p and -Eaip 2 , at the associated stretch directions are 
given by the identities: 

_^lL + ^ii = 2 , ^11L + = 2. (5.2) 

Eu Pi v\ Eiip 2 V2 — 1 

The first identity follows from the pair of equations (3.13) and the second from 
(3.17). Equations (5.2) indicate that if v\ or vi become large in magnitude then 
the second term in the left member is negligible, and the associated value of the 
Young's modulus is approximately one half of the value in the 111 direction. Thus, 
large values of \v\ occur in directions at which E sa |£ nl . Such directions, by their 
nature, are close to 111. 




-1.5 -1 -0.5 0.5 1 1.5 2 2.5 

001 



Figure 5. The global maximum of Poisson's ratio based on table 2. The value of v max 
depends upon the location of (fooi, ^no) in five distinct regions defined by the heavy 
lines. The dashed line delineates the (small) region in which tWx > 2, from eq. (5.1). 

We note that v\ appears in both figures 4 and 5. The double occurrence is not 
surprising if one considers that fooi, ^110 and i>o also occur in both the minimum 
and maximum. It can be checked that in the region where v\ is the maximum value 
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in figure 5 it satisfies —1 < v\ < —1/2. In fact it is very close to but not equal to 
fil in this region, and numerical results indicate that — Viiq\ < 4 x 10~ 4 in 
this small sector. 

What is special about the transition values in figures 4 and 5: v 00 i = 3/2 and 
Vxiq = —1/2? Quite simply, they are the values of V\ and v 2 as the stationary 
directions n approach the face diagonal direction 110. Thus, V\ and v 2 are both 
the continuation of the face diagonal value ^j , but on two different branches. See 
Appendix B for further discussion. 




-1 -0.5 o V 0.5 1 1.5 2 

001 



Figure 6. The 44 materials considered are indicated by dots on the chart showing the 

Vmi n regions, see figure 4. 



(b) Application to cubic materials 

We conclude by considering elasticity data for 44 materials with cubic symmetry, 
figure 6. The data are from Musgrave (2003) unless otherwise noted. The cubic 
materials in the region where v m - m = fooi are as follows, with the coordinates 
Oooi,fiio) for each: GeTeSnTef (mol% GeTe=0) (0.01, 0.70), RbBrf (0.06, 0.64), 
KI (0.06, 0.61), KBr (0.07, 0.59), KC1 (0.07, 0.56), Nbf (0.21, 0.61), AgCl (0.23, 
0.61), KF1 (0.12, 0.49), CsCl (0.14, 0.44), AgBr (0.26, 0.55), CsBr (0.16, 0.40), NaBr 
(0.15, 0.38), Nal (0.15, 0.38), NaCl (0.16, 0.37), CrVf (Cr0.67 at.% V) (0.15, 0.35), 
Csl (0.18, 0.38), NaFl (0.17, 0.32). This lists them roughly in the order from top 
left to lower right. Note that all the materials considered have positive i'ooi- The 
materials with i/ m ; n = f (n a l so have ^ max = ^iio 7 so the coordinates of the above 
materials correspond to their extreme values of v. The extreme values are also given 
by the coordinates in the region with v^ lVL = Vi%q, i/ max = vooi- The materials there 

f Data from Landolt and Bornstein (1992), see also Cazzani and Rovati (2003). 
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Table 3. Properties of the 11 materials of cubic symmetry in figure 6 with u x \ a < —1/2. 
The boldfaced numbers indicate i/ min and f m „. Unless otherwise noted the data are from 
Landolt and Bornstein (1992). G&S indicates Gunton and Saunders (1975). 



Material 


^001 


^110 




Pi 


V2 


P2 


Ml//"2 


/3-brass (Musgrave 2003) 


1.29 


-0.52 


-0.52 


0.15 






8.5 


Li 


1.29 


-0.53 


-0.54 


0.21 






8.8 


AINi (at 63.2% Ni and at 273 K) 


1.28 


-0.55 


-0.55 


0.25 






9.1 


CuAINi (Cul4% A14.1 % Ni) 


1.37 


-0.58 


-0.59 


0.32 






10.2 


CuAINi (Cul4.5% A13.15% Ni) 


1.41 


-0.63 


-0.66 


0.42 






12.1 


CuAINi 


1.47 


-0.65 


-0.69 


0.45 






13.1 


AINi (at 60% Ni and at 273 K) 


1.53 


-0.68 


-0.74 


0.50 


1.53 


0.18 


15.0 


InTl (at 27% Tl, 290K) (G&S) 


1.75 


-0.78 


-0.98 


0.62 


1.89 


0.59 


24.0 


InTl (at 28.13% Tl) 


1.78 


-0.81 


-1.08 


0.66 


2.00 


0.63 


28.6 


InTl (at 25% Tl) 


1.82 


-0.84 


-1.21 


0.70 


2.14 


0.68 


34.5 


InTl (at 27% Tl, 200K) (G&S) 


1.93 


-0.94 


-2.10 


0.83 


3.01 


0.82 


90.9 



are: Al (0.41, 0.27), diamond (0.12, 0.01), Si (0.36, 0.06), Ge (0.37, 0.02), GaSb 
(0.44, 0.03), InSb (0.53, 0.03), CuAuf (0.73, 0.09), Fe (0.63, -0.06), Ni (0.64, -0.07), 
Au (0.88, -0.03), Ag (0.82, -0.09), Cu (0.82, -0.14), a-brass (0.90, -0.21), Pbf (1.02, 
-0.20), Rbf (1.15, -0.40), Csf (1.22, -0.46). 

Materials with < —1/2 are listed in table 3. These all lie within the region 
where the minimum is u\, and of these, five materials are in the sub-region where 
the maximum is v^. Three materials are in the sub- regions with v\ < — 1 and 
V2 > 2. These Indium Thallium alloys of different composition and at different 
temperatures are close to the stability limit where they undergo a martensitic phase 
transition from face-centered cubic form to face-centered tetragonal. The transition 
is discussed by, for instance, Gunton and Saunders (1975), who also provide data 
on another even more auxetic sample: InTl (at 27% Tl, 125K). This material is so 
close to the ^2 = vertex, with j/ oi = 1.991, v x \^ — —0.997 and ^1/112 = 1905 (!) 
that we do not include it in the table or the figure for being too close to the phase 
transition, or cquivalently, too unstable (it has v\ = —7.92 and V2 = 8.21). 

We note that the stretch directions for the extremal values of v, defined by 
n = llpi and n = llp2, are distinct. As the materials approach the ^2 = vertex 
the directions coalesce as they tend towards the cube diagonal 111. The three 
materials in table 3 with ^ m ; n < — 1 and ^ max > 2 are close to the incompressibility 
limit, the line k = oo in figure 3. In this limit both the cube diagonal and axial 
Poisson's ratios tend to 1/2, i.e. vm = vo = 1/2, and 



1 1 

4~ 1" 



V2 



4' 



Pi =Vi 



'1-3, 



k — ► oo . (5.3) 



These are reasonable approximations for the last three materials in table 3, which 
clearly satisfy v\ + V2 ~ 1 , and pi wp 2 . 



Article submitted to Royal Society 



Poisson's ratio in cubic materials 



17 



6. Summary 

Figures 4 and 5 along with tables 1 and 2 are the central results which summarize 
the extreme values of Poisson's ratio for all possible values of the elastic parameters 
for solids with positive strain energy and cubic material symmetry. The application 
of the related formulas to the materials in figure 6 shows that values less than — 1 
and greater than +2 are associated with certain stretch directions in some Indium 
Thallium alloys. 

Discussions with Prof. T. C. T. Ting are appreciated. 



Appendix A. Extreme values of D(n, m) for a given n 

The extreme values of Z?(n, m) as a function of m for a given direction n can be 
determined using Lagrange multipliers A, p, and the generalized function 



/(m) = D(n, m) — A|m| 2 — 2pn ■ m . 



(Al) 



Setting to zero the partial derivatives of / with respect to m\, mz, m.3, implies 
three equations, which may be solved to give 



pn 2 



pn 3 



v n 2 - A' n\ — A' n\ - X 1 



(A 2) 



where A, p follow from the constraints n • m = and |m| 2 = 1. These are, respec- 
tively, (2.16) and 



(n 2 -X) 2 (n 2 -A) 2 (n 2 -X) 2 



1. 



(A3) 



Equation (2.16) implies that A is a root of the quadratic equation (2.17) and (A3) 
yields the normalization factor p. These results are summarized in equations (2.18) 
and (2.19). 

It may be easily checked that the generalized function / is zero at the extremal 
values of D. But / = D — A, and hence the extreme values of D(n, m) are simply 
the two roots of the quadratic (2.17), < A_ < A + < 1/2. Note that the extreme 
values depend only upon the invariants of the tensor M with components My = 
Dijkirikni. Although this is a second order tensor and normally possesses three 
independent invariants, one is trivially a constant: trM = 1. The others are, e.g. 
trM 2 = n\ + n\ + n i = 1 ~ 2F(n) (see equation (2.8)) and detM = n\n\n\. 

The above formulation is valid as long as {n\ — n 2 .)^ 2 . — n§)(n§ — n\) ^ 0. For 



v 2 , then A_, Aj 



2A"2 

min, max (n 2 , 3n 2 n 2 ). The m vector associated 

2 _^ ^2 



instance, if 

with A = n\ is undefined, according to (A 2). However, by taking the limit n'2 — >• n\ 
it can be shown that m — > ±(1, — 1, 0)/y/2. The other vector corresponding to 
A = 3n 2 n 2 has no such singularity, and is m = ±(713, n.3, — 2ni)/v / 2- 

The identity (3.4) may be obtained by noting that each term can be split, e.g. 
n i/( n i — A) = n 2 + X/(n 2 — A), then using the fundamental relation (2.16) with 
n\ + n 2 + n\ = 1. Various other identities can be found, e.g. 



A+)(n 2 -A_) 



A+)K-A_) 



(nl 



A+)(n?-A_) 



(A 4) 
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Appendix B. Analysis 

Here we derive stationary conditions for directions n along the edges of the triangle 
in figure 2 by direct analysis. Numerical tests are performed for the entire range of 
material parameters. The results are consistent with and reinforce those of §3. 

The limiting Poisson's ratios of (4.4) are expressed f±(n) = v±{ot,(3) in terms 
of two numbers, where 

n? = (l + a)(l-/3), nI = (l-a)(l-/3), n§ = \ + 2/3. (B 1) 

The range of (a, (3) which needs to be considered is < a < 1, 3 ^ 3 ° a ^ < (3 < 
1/3, corresponding to the triangle in figure 2. This parameterization allows quick 
numerical searching for global extreme values of v for a given cubic material. 

We first consider the three edges as shown in figure 2 in turn. Edge 1 is defined 
by a = 1, 1/12 < < 1/3. The limiting values are V-(l,f3) = vq/[\-(\ - 0){\ + 
2/3)2xo] and ^+(l,/3) = 1 — z^_(l,/3)(l — ^o)/^o- The extreme values are obtained 
at the ends: f-(l, 1/12) = ^ooi, ^+(1, 1/12) = vno, 1/3) = 1/3) = v^. 

These possible global extreme values agree with those of §3. 

Inspection of figure 2 shows that edges 2 and 3 can be considered by looking at 
v±(0, (3) for — 1/6 < (3 < 1/3. Straightforward calculation gives 

v-{y,p)-z — 7i — — , — - — — . [a 2) 

l-(^-3/3 2 )xo l-(i-3/3 2 )xo 

A function of the form f /g is stationary at f/g = fig'. Applying this to the ex- 
pressions in (B 2) implies that the extreme values of v_ and v + satisfy, respectively, 

"-(0,/3) = ^, " + M = 1~. (B3) 

Combining equations (B 2) and (B 3) gives in each case a quadratic equation in (3. 
Thus, the extreme values of v_ and v + are at (3 = (3_± and (3 — (3+±, the roots 
of the quadratic equations. The first identity, (B3)i was found by Ting & Chen 
(2005), their equation (4.15). 

To summarize the analysis for the three edges: Extreme values of Poisson's 
ratio on the 3 edges are at the ends of edge 1, and on edges 2 and 3 given by v± of 
equations (B2)-(B3). 



(a) Numerical proof of tables 1 and 2 

A numerical test was performed over the range of possible materials. This re- 
quired searching the entire two-dimensional range for a, (3. Consideration of all 
possible materials then follows by allowing the material point to range throughout 
the triangle of figure 3. In every case it is found that the extreme values of v occur 
on the edge of the irreducible l/48th element of the cube surface. Furthermore, the 
extreme values are never found to occur along edge 2. Extreme values on edge 3 in 
figure 2 can be found by considering edge 3' instead, i.e. a = 0, — 1/6 < (3 < 0. This 
implies as possible extrema one of V-(0,(3-±) and one of i/+(0, /3+±). We define 
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these as v[ = z/_(0, /3 ), and v' 2 = i/ + (0,/3 H ), where the signs correspond to the 

sign of the discriminant in the roots, then they are given explicitly as 




It may be checked that u[ = v\ and v' 2 = v 2 , in agreement with equation (3.19). 

The numerical results indicate the potential extrema come from the five values: 
^o, ^ooi) u i an d v i- K turns out that each is an extreme for some range of 

material properties. Thus, the first four are necessary to define the global minimum, 
see table 1 and figure 4, while all five occur in the description of the global maximum, 
in table 2 and figure 5. 

Although a mathematical proof has not been provided for the veracity of tables 
1 and 2, and figures 4 and 5, it is relatively simple to do a numerical test, a poste- 
riori. By performing the numerical search as described above, and subtracting the 
extreme values of tables 1 and 2, one finds zero, or its numerical approximant for 
all points in the interior of the triangle of possible materials, figure 2. 

In order to further justify the results as presented, the next subsection gives 
arguments for the occurrence of the special values —1/2 and 3/2 in figures 4 and 5. 



(b) Significance of — ^ and | 

Suppose Poisson's ratio is the same for two different pairs of directions: v(n, m) = 
i/(n*,m*). The pairs (n, m) and (n*,m*) must satisfy, using (4.3), 

D(n,m)-D(n'm') = 

F(n)-F(n*) v ' 

For instance, let n* = 110, m* = 001, so that v — ^ooi- Equation (B6) implies 
that the same Poisson's ratio is achieved for directions (n, m) satisfying 

^ (n ' m) -.oor. (B7) 



2F(n) - 1/2 



Note that this is independent of v and xo- We choose Vqoi specifically because it has 
been viewed as the candidate for largest Poisson's ratio, until Ting & Chen (2005). If 
it is not the largest, then there must be pairs (n, m) other than (110, 001) for which 
(B 7) holds. However, it may be shown using results from §2 that the minimum of 
the left member in (B7) is 3/2, and the minimum occurs at n = 110, as one might 
expect. This indicates that ^ooi must exceed 3/2 in order for the largest Poisson's 
ratio to occur for n other than the face diagonal 110. 

Returning to (B6), let v = vn a , then i/(n,m) = vn if 



i - D(n, m 



2F(n) 



=-i/il . (B8) 



2 



Using equation (2.8) and the previous result, it can be shown that the minimum of 
the left member in (B8) is 1/2, and the minimum is at n = 110. Hence, Vii must 
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be less than —1/2 in order for the smallest Poisson's ratio to occur for n other than 
the face diagonal 110. These two results explain why the particular values ^ooi = § 
and = — i appear in tables 1 and 2 and in figures 4 and 5. 
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